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theory free of the sign problem such as two-color QCD at finite chemical potential. 
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1. Introduction 

The investigation the phase diagram of QCD in the temperature-chemical potential plane has 
a deep relevance and implications on cosmology, astrophysics and in the phenomenology of heavy 
ions collisions. The lattice formulation of QCD is the only tool to approach this important issue 
starting from first principles. Unfortunately the study of QCD at non-zero baryonic density by nu- 
merical simulations on a space-time lattice is plagued by the well-known sign problem: the fermion 
determinant is complex and the Monte Carlo sampling becomes unfeasible. One of the possibilities 
to circumvent this problem is to perform Monte Carlo numerical simulations for imaginary values 
of the baryonic chemical potential, where the fermion determinant is real and the sign problem is 
absent, and to infer the behavior at real chemical potential by analytic continuation. The method 
of analytic continuation [jl|-[l^] is well-founded and works fine within the limitations posed by the 
presence of non-analyticities and by the periodicity of the theory with imaginary chemical poten- 
tial [0]. 

It is very important to answer the question about which is the optimal way to extract infor- 
mation from data taken at imaginary values of the chemical potential. This is equivalent to an- 
swer which is the best interpolating function for data at imaginary chemical potential that could 
analytically continued in order to get physical predictions for real values of fi. The aim of our 
investigations in the past two years has been to to study limitations and possible improvements of 
the method of analytic continuation []|, [jj, 19]. In order to study this problem we have considered 



SU(2) (two-color QCD) and SU(3) at finite isospin. Indeed these theories are free of the sign prob- 
lem and Monte Carlo numerical simulations at real values of the chemical potential or at real values 
of isospin potential are feasible. Therefore it is possible to compare the analytic continuations with 
the data from direct simulations allowing at the same time to discriminate between interpolating 
functions an to test the range of the reliability of the analytic continuation. Here we briefly review 
results obtained in studying analytic continuation of physical observables and of the critical line in 



two-color QCD. Results obtained for SU(3) at finite isospin has been reviewed in ref. [2C]. 
2. Analytic continuation of physical observables 



As already shown long ago Q17[], the partition function of any SU(N) gauge theory with non- 
zero temperature and imaginary chemical potential, /I = is periodic in 6 = \ijjT with period 
2n/N and that the free energy F is a regular function of 8 for T < Te, while it is discontinuous at 
6 = 2n(k + 1 /2)/N, k = 0, 1,2, . . ., for T > Te, where Te is a characteristic temperature, depending 
on the theory. 

We have considered SU(2) in presence of tif = 8 degenerate staggered fermions of mass 
am = 0.07. Figure 1 shows the the tentative phase diagram in the (/X/,jS) plane for this theory 
(in correspondence of a fermion mass am = 0.07) with Pe — 1-55 [ [HI |22| ] and j8 c ~ 1.41 [p3|]. We 
have performed numerical simulations on a 16 3 x 4 lattice using the exact <p algorithm [ pT| ] with 
dt = 0.01 (typical statistics: 20k trajectories). We have performed a careful test of the analytic 
continuation of physical observables. A detailed discussion of the results obtained is reported in 
ref. [Q]. In order to show the importance of a careful choice of the interpolating function of the 
imaginary chemical potential data, in figure ^| we display results for two different observables in 
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Figure 1: Phase diagram in the (/3 , /i/)-plane; N is the number of colors, N r the extension of the lattice in the 
temporal direction. The numerical values for $e and j3 c are valid for SU(2) in presence of nf = 8 degenerate 
staggered fermions with mass am = 0.07 



correspondence of two different values of j8. We have shown that the use of ratio of polynomials as 
interpolating function can lead to a dramatic improvement in the analytic continuation of physical 
observables. 

The aim of our subsequent investigations has been to understand if a careful choice of the 
interpolating function can also improve the continuation of the critical line. 
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Figure 2: (Left) Negative side of the horizontal axis: imaginary part of the fermion number density vs. the 
imaginary chemical potential at j3 = 1.45. Positive side of the horizontal axis: real part of the fermion num- 
ber density vs. the real chemical potential at j3 = 1 .45. The green (blue) solid lines represent the polynomial 
(ratio of polynomials) interpolating function; the dashed lines give the corresponding uncertainty, coming 
from the errors in the parameters of the fit. (Right) Chiral condensate vs. ju 2 at /3 = 1.90. The green (blue) 
solid lines represent the polynomial (ratio of polynomials) interpolating function; the dashed lines give the 
corresponding uncertainty, coming from the errors in the parameters of the fit. 



3. Analytic continuation of the critical line 

The determination of the critical line in the (T,fi) plane is of overwhelming importance for the 
study of strong interactions at finite temperature and baryon density. The analytic continuation of 
the (pseudo-)critical line on the temperature-chemical potential plane is well-justified, but a careful 
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test in two-color QCD and three-color QCD with finite isospin chemical potential has cast some 
doubts on its reliability [|l8|-[20|]. 
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Table 1: Summary of the values of j8 c (fi 2 ) obtained by fitting the peaks of the susceptibilities of chiral 
condensate (Tpy), Polyakov loop (L) and plaquette (P) in SU(2) on a 16 3 x 4 lattice with fermion mass 
ara=0.07 ', For each interpolation the ^ 2 /d.o.f. is given. 

In this section we present our results in the determination of the critical line in two-color QCD 
using the method of analytic continuation. Contrary to the case of physical observables discussed 
in the previous section, the theoretical basis is not straightforward since it relies on the assumption 
that susceptibilities, whose peak signals the presence of the transition, be analytic functions of the 
parameters on a finite volume ^ We have tested the method of analytic continuation in the case 
of two-color QCD and in the case of QCD ad finite isospin density [18, As for usual QCD 



simulations, we have determined the critical line for imaginary values of the chemical potential 
(/J. 2 < 0) and interpolate them by suitable functions to be continued to pr > 0. In order to test the 
reliability of the analytic continuation, the prediction obtained at real /J. has been compared with 
direct determinations of the transition line. 
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Figure 4: Critical couplings obtained from the susceptibility of chiral condensate in SU(2) on a 16 3 x4 
lattice with ara=0.07 (left) and am=Q.2 (right), together with a linear fit (dotted line) in (ajj.) 2 to data with 
ji 2 < 0. The solid lines around the best fit line delimit the 95% CL region. 
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Figure 5: Critical couplings obtained from the susceptibility of Polyakov loop in SU(2) on a 16 3 x 4 lattice 
with am=0.07 (left) and am=0.2 (right), together with a linear fit (dotted line) in (ajx) 2 to data with [I 2 < 0. 
The solid lines around the best fit line delimit the 95% CL region. 



It should be remarked that on a finite volume there are no true nonanalyticities, therefore the 
location of the critical line may depend on the observable chosen to probe the transition. Conse- 
quently we have determined the (pseudo-)critical coupling j8 c (/i 2 ) by looking at the peaks of the 
susceptibilities of three different observables: the chiral condensate, the Polyakov loop, and the 
plaquette. In figure || we show an example of this determination. We have fitted the peak according 
to a Lorentzian function. The result obtained in this way for j6 c (jU 2 ) agrees well with the result 
obtained by means of the multihistogram reweighting. 

In table [T] the values of j3 c (;lt 2 ) obtained in correspondence of the different "probe" observable 



5 



Analytic continuation in QCD-like theories 



L. Cosmai 



(at fermion mass am = 0.07) are shown. They depend very weakly on the observable considered 
(only in one case the relative deviation between two determinations at the same \i 2 slightly ex- 
ceed 3a). The strategy is now to interpolate the critical j8's obtained at fixed imaginary chemical 
potential with an analytic function of /x, to be then extrapolated to real chemical potential. For a 
theory free from the sign problem (such as two-color QCD) the extrapolated curve can be com- 
pared with the determinations of the critical j8 's at real chemical potential obtained through a direct 
computation. 

As displayed on the left sides of figure |] and figure |5[ imaginary chemical potential data 
are very well fitted with a linear polynomial in /I 2 . If different functions are used (larger order 
polynomials, ratio of polynomials) the fit puts to value compatible with zero all parameters (except 
two of them) thus reducing again to a first order polynomial in jU 2 . This is in marked difference 
with what we found (see previous section and ref. [Q]) for the analytic continuation of physical 
observables, where the ratio of polynomials performed very well. Moreover, on the left sides of 
figures || and ||, we can clearly see a quite significative deviation between extrapolation and direct 
determination of the critical line at real chemical potential. The discrepancy found could imply 
that or the the critical line is not analytic on the whole interval of 2 or that the interpolation at 
/I 2 < is not accurate enough to correctly reproduce the behavior at ji 2 > 0. Indeed we have found 
(see figure ||) that a polynomial of third order in jj. 2 nicely fits all data for j6 c (jU 2 ) and therefore 
the critical line is analytic on the whole interval of \i 2 . A possible conclusion is that for /I 2 < 
the jit 4 and /x 6 terms compensate each other at large negative values of /I 2 so that the effective 
interpolating function of the data at p} < is a first order polynomial in /I 2 while for /I 2 > the 
/I 4 and jit 6 terms work in the same direction and their contribution cannot be neglected. 

In order to verify if this scenario is peculiar to SU(2) we have investigated the same theory 
with a different mass am = 0.2 for the «/ = 8 degenerate fermions. We have also examined the 
case of SU(3) at finite isospin [19, 2(J. The reason for consider SU(2) at a larger value for the quark 
mass is that at any fixed temperature T the critical value of /I gets larger, consequently the critical 
line could become less curved in the physical region /I 2 > 0. Accordingly higher order terms in \i 2 
in the description of the critical line by a polynomial could be less important. We have sampled the 
critical line for SU(2) and quark mass am = 0.2. The best interpolation of /3 c (jU 2 ) data at jj, 2 < is 
a polynomial linear in \i 2 and, at variance with the case of quark mass am = 0.07, the extrapolation 
to /I 2 > compares very well with the direct determination of j3 c (^t 2 ) in that region. So we can 
argue that the the extrapolation to \i 2 > works definitely better for larger quark masses, i.e. away 
form the chiral limit. 



4. Conclusions 

We have reported here an investigation which is part of a larger project devoted to study the 
method of analytic continuation in QCD-like theories free from the sign problem and to improve 
its predictivity in view of its application to QCD. 

For what concerns analytic continuation of physical observables we have shown that a con- 
siderable improvement can be achieved, when extrapolating data from imaginary to real chemical 
potentials, if ratios of polynomials are used as interpolating functions (for a thorough discussion 
see ref. [§]). 
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Figure 6: Same as in figure 4, but with results of a fit to all data including term up to order /I 



We have also presented results for the analytic continuation of the critical line in the (r,/i) 
plane from imaginary to real chemical potential both in the case of two-color QCD. We have found 
that the critical line around pL = can be described by an analytic function. Indeed, a third order 
polynomial in /I 2 fits all the available data for the critical coupling. 

We have shown that there is a clear indication that in the chiral limit high-order terms in the 
polynomial interpolation play a relevant role at \i 2 > but are less visible at /I 2 < 0, this calling 
for extremely high accuracy in detecting such terms from simulations at /I 2 < 0. The predictions 
for the pseudocritical couplings at real chemical potentials may be wrong if data at imaginary /I 
are fitted according to a linear dependence. 

All the issues above have undergone further investigation in a different theory such as SU(3) 
at finite isospin (results for SU(3) at finite isospin are discussed in refs. [19|, 2(J). 

The lessons we learned in studying analytic continuation in QCD-like theory free from the 
sign problem will be used in the near future to the determination the critical line for real QCD. 
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